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1.  Background 

The  general  purpose  of  this  research  is  to  study  the  physics  and  dynamics  of 
stratiform  rain.  Our  contribution  to  the  project  was  to  acquire  vertical  profiles  of 
Doppler  velocity  and  signal  intensity  with  a  94  GHz  (3  mm  wavelength)  Doppler 
radar  operated  in  a  vertically  pointing  mode  and  to  analyze  the  94  GHz  radar 
observations  together  with  joint  Ka  radar  observations  and  rain  measurement  at 
the  ground  using  a  raingage  and  a  disdrometer.  Using  a  94  GHz  Doppler  radar 
for  such  a  project  yields  a  very  high  sensitivity  on  weak  targets,  a  high  spatial 
resolution  due  to  the  availability  of  a  narrow  beam  and  a  short  pulse  width. 
Also,  the  raindrops’  bacl;sca‘tering  cross  section  vs  size  relaiionsmp  at  94  GHz 
exhibits  deep  oscillations  within  the  raindrop  size  range  which  provides  a  means 
to  identify  sizes  in  a  Doppler  spectrum  observed  at  vertical  incidence. 

Our  effort  was  divided  into  the  following  parts: 

•  Preparation  of  the  94  GHz  radar 

•  Pilot  experiment  in  May-June  1987 

•  Second  field  experiment  in  November-December  1987 

•  Data  analysis  and  report  preparation 


2.  94  GHz  Doppler  radar 

The  radar  and  its  performance  during  various  projects  was  presented  else¬ 
where  [Lhermitte,  1987,  1988a].  A  photograph  of  the  radar  is  shown  in  Fig.  1 
and  the  main  characteristics  of  the  radar  when  used  in  the  GL  project  are  the 
following: 

•  Peak  power:  1.2  kW 

•  Receiver  noise  figure:  6.5  db  DSB 

•  Minimum  detectable  signal:  -99  dBm  (-115  dBm  with  3  s  integration) 
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•  Two  (transmit  and  receive)  3-foot  antennae. 

The  radar  reflectivity  of  precipitation  or  clouds,  g,  can  be  expressed,  knowing 
the  intensity,  Pr,  of  the  backscattering  signal  they  produce,  by  the  expression: 


10  log  g  =  10 log  PT  —  10 log  Pt  +  lOlog  4 nR2  —  10 log  h  —  10 log  Ae  (1) 

where  R  is  the  target  distance,  Pt  is  the  transmitter  peak  power,  h  is  the  pulse 
volume  radial  dimension,  and  Ae  the  antenna  effective  area. 

With  the  radar  characteristics  above  and  including  a  0.6  (two-way)  antenna 
efficiency  coefficient,  rj  is  given  by: 

10  log  x]  =  PT(dBm)  -  23.5  +  20 log  R  (2) 

where  g  is  in  cm-1  and  R  is  in  kilometers.  In  the  case  of  a  Rayleigh  target  (hy¬ 
drometeors  having  a  size  much  smaller  than  the  radar  wavelength  and  practically 
applicable  only  to  cloud  droplets),  we  have: 

dBZ  =  Pr(dBm)  +56  +20  logR  (3) 

where  dBZ  ~  lOlogZ  is  the  reflectivity  factor  expressed  in  units  familiar  to  radar 
meteorologists. 

The  radar  is  equipped  with  a  processor  delivering  profiles  of  mean  Doppler 
and  mean  signal  intensity  estimates  simultaneously  at  200  range  gates  spaced 
by  400  ns  (60  m).  For  the  November  1987  experiments,  a  new  recording  unit, 
capable  of  storing  4096  2-byte  samples  of  the  received  signal  (I  and  Q)  at  a 
selectable  range  gate,  was  also  included  in  the  radar  equipment  so  that  complete 
Doppler  spectra  could  be  calculated  at  any  of  the  range  gate  positions.  The  data 
are  acquired  successively  at  16  range  gates  as  the  range  is  stepped  by  200  ns  (30 
m)  intervals,  thereby  allowing  observation  of  480  m  sections  of  Doppler  spectra 
vertical  profiles. 

3.  Millimeter  wave  scattering  and  absorption 

At  94  GHz,  scattering  and  absorption  of  millimeter  wave  radiation  by  rain¬ 
drops  is  determined  using  Mie  scattering  functions  rather  than  the  more  familiar 
(and  simpler)  Rayleigh  scattering  applicable  to  centimeter  wave  radiation. 
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This  is  illustrated  by  Fig.  2  which  shows  the  backscattering  cross  section 
of  water  spheres,  at  both  35  GHz  and  94  GHz,  as  a  function  of  their  diameter, 
D.  The  Mie  function  at  94  GHz  exhibits  a  maximum  at  approximately  1.1  mm 
diameter  and  a  well  defined  minimum  at  1.57  mm.  Fig.3  illustrates  the  difference 
between  backscattering  cross  sections  for  ice  and  water  at  the  94  GHz  frequency. 
Note  that,  while  liquid  water  drops  are  stronger  scatterers  than  ice  spheres  within 
the  Rayleigh  assumption,  this  does  not  hold  for  particles  not  small  with  respect 
to  the  radar  wavelength,  since,  in  that  case,  ice  spheres  are  better  scatterers  than 
water  drops  of  the  same  size. 

Signal  attenuation  is  also  much  more  significant  at  millimeier  waves,  with 
both  clouds  and  precipitation  being  strong  absorbers.  Within  the  Rayleigh  scat¬ 
tering  assumption  (applicable  to  cloud  droplets  even  at  94  GHz),  signal  absorp¬ 
tion  is  proportional  to  liquid  water  and  is  controlled  by  the  imaginary  term  of 
the  water  index  of  refraction  which  varies  with  droplet’s  temperature.  This  is 
illustrated  by  Fig.  4  which  shows  the  94  GHz  absorption  coefficient  of  cloud  or 
fog  (per  g  m~ 3)  as  a  function  of  water  temperature  [Lhermitte,  1986]. 

At  94  GHz,  the  Rayleigh  assumption  is  not  applicable  to  precipitation  parti¬ 
cles  and  signal  attenuation  by  precipitation  must  be  evaluated  using  the  Mie  func¬ 
tions.  Theoretical  attenuation  coefficients  computed  using  this  method  [Lher¬ 
mitte,  1990]  are  shown  in  Fig.  5  together  with  an  absorption  coefficient  vs  rain 
intensity  relationship  proposed  by  Wallace  [19S8],  as  a  least  squares  fit  to  his 
experimental  data.  Also  shown  in  Fig. 5  is  a  measurement  (circled  cross)  based 
on  experimental  data  acquired  and  interpreted  as  part  of  the  research  performed 
within  the  GL  contract.  These  experimental  data  indicate  that  the  attenuation 
calculated  from  extinction  Mie  cross  sections  is  greater  than  that  actually  shown 
by  measurements.  This  may  be  due  to  the  contribution  of  forward  scattering  to 
the  radiation  propagating  forward  in  the  scattering  medium.  Investigations  of 
the  magnitude  of  this  effect  was  beyond  the  scope  of  the  present  effort.  However 
such  investigations  are  possible  using  precipitation  models  and  the  Mic  fun 
tions.  They  are  recommended  to  provide  the  theoretical  background  needed  to 
fully  understand  signal  attenuation  at  millimeter  waves. 

4.  Pilot  experiment 

Two  field  experiments  (May- June  and  November-December  1987)  were  per¬ 
formed  at  the  GL  site  in  Sudbury,  Mass.  The  first  experiment  was  devoted  to 
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building  the  radar  shelter,  testing  the  radar,  and  acquiring  observations  primar¬ 
ily  for  the  purpose  of  evaluating  the  radar  performance  and  the  signal  and  data 
processing  methods.  Some  observations  in  May  were  conducted  with  the  radar 
moved  outside  and  covered  with  a  waterproof  protective  canvas. 

The  shelter  was  designed  for  housing  of  the  radar  and  associated  equipments 
in  all  weather  conditions  and  was  planned  to  be  installed  at  a  selected  location 
near  the  main  GL  Sudbury  facility  in  November.  It  was  made  of  lumber  and 
plywood  and  covered  with  corrugated  fiberglass  panels.  The  roof  material  was 
tested  for  94  GHz  rodiatior  attenuation  and  it  was  found  that,  if  dry,  this  material 
produces  a  2-3  dB  two-wa>  '  tenuation  which  was  judged  acceptable.  The  actual 
use  of  the  shelter  in  November  showed  that  this  building  provided  an  adequate, 
waterproof,  shelter  in  snow  or  rain  conditions.  However,  in  these  conditions,  the 
roof  would  hold  water  or  snowflakes  which  would  drastically  increase  the  signal 
attenuation  produced  by  that  roof. 

Suitable  stratiform  rain  conditions  occurred  on  May  15.  A  well  defined  melt¬ 
ing  level  indicated  by  a  sudden  and  drastic  increase  of  mean  Doppler  due  to 
melting  of  snowflakes  as  they  fall  through  the  0°C  isotherm  was  observed.  The 
presence  of  a  well-  defined  melting  level  suggests  the  absence  of  significant  air 
vertical  velocity,  but  we  found  that,  in  the  rain  below  the  melting  level,  some 
convection  (50  cm  s-1  up-downdrafts  indicated  by  time  variability  of  the  mean 
Doppler  velocity  discussed  below t  may  have  existed  at  certain  times. 

The  data  recording  started  at  09:3G  EDT.  The  radar  data  indicated  a  cloud 
top  at  approximately  8  km  but  there  was  no  precipitation  observed  at  the  ground. 
At  09:45  very  sparse  drops  fell  but  the  rain  was  very  light  and  not  measurable. 
The  precipitation  then  increased  steadily  to  reach  drizzle  to  light  rain  (approxi¬ 
mately  1  mm  /fir-1)  conditions  later.  The  observations  lasted  slightly  over  one 
hour  during  which  more  than  900  vertical  profiles  of  mean  Doppler  and  signal 
intensity  were  obtained. 

A  typical  example  of  the  profiles  of  mean  Doppler  and  signal  intensity  is 
shown  in  Fig.  6.  The  layer  separating  dry  snowflakes  above  and  raindrops  below 
is  clearly  shown  in  the  mean  Doppler  (and  to  a  lesser  degree  signal  intensity)  ver¬ 
tical  profiles.  Four  successive  profiles  of  both  signal  intensity  and  mean  Doppler, 
taken  at  approximately  5s  time  intervals,  are  shown  in  Fig.  7.  The  signal  inten¬ 
sity  shown  in  Figs.  6  and  7  (and  Fig.  8  below)  is  corrected  for  f?2  and  signal 
attenuation  estimated  from  the  Wallace  [1988]  data  and  the  rain  intensity  at 
the  ground.  Also,  the  velocity  profiles  are  corrected  for  air  density,  p,  using  the 
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relationship  V(D)  =  V0(D)  (p/ p0)x  [Foote  and  DuToit,1969j,  where  V0 (D)  is  the 
drops’  terminal  velocity  at  the  ground  and  x  =  0.45  (see  equation(4)  below). 
Fig. 7  reveals  a  remarkable  steadiness  of  the  vertical  profiles  (especially  signal 
intensity). 

The  most  striking  feature  in  the  data  presented  in  Figs.  6  and  7  is  indeed 
the  drastic  variation  with  height  of  the  mean  Doppler  velocity  produced  by  the 
melting  of  snowflakes  when  passing  through  the  melting  layer.  There  is  a  large 
increase  (1.5  m/s  to  4.5  m/s)  of  the  mean  Doppler  in  a  200  m  (approximately 
2.9  km  to  2.7  km)  vertical  interval.  This  reflects  the  acceleration  of  precipitation 
particle  fall  speed  as  they  evolve  from  an  ice-crystal-snowflake  nature  to  a  rain¬ 
drop  shape.  This  is  similar  to  what  was  observed  much  earlier  with  centimeter 
wave  meteorological  radars. 

However,  centimeter  wave  radars  observations  revealed  a  well-  defined  max¬ 
imum  of  radar  reflectivity  (from  which  the  name,  ’bright  band’,  was  derived) 
occurring  just  below  the  0 °C  level  [Lhermitte  and  Atlas,  1963].  This  wras  ex¬ 
plained  earlier  by  a  combination  of:  i. )  above  the  maximum,  an  increase  of  radar 
reflectivity  due  to  the  change  of  index  of  refraction  from  ice  to  w'ater  (7-S  dB  at 
centimeter  wavelengths);  ii)  below  the  maximum,  a  decrease  of  the  concentra¬ 
tion  of  precipitation  particles  due  to  their  acceleration  during  melting  and  the 
reduced  particle  size  after  melting.  Analysis  of  ’bright  band’  data  were  based  on 
the  fact  that,  a  centimeter  wavelengths  the  precipitation  particles'  radar  cross 
section  is  evaluated  using  the  Rayleigh  assumption,  e.  g.  the  radar  cross  sec¬ 
tion  always  increases  monotonically  with  particle  size.  Also,  it  was  assumed  that, 
each  particle  retains  its  identity  wdth  no  mutual  interactions,  break-up  or  growth 
when  it  evolves  from  an  ice  crystal-snowflake  shape  to  its  final  raindrop  shape. 
However,  aggregation  and  break-up  may  occur  and  can  appreciably  modify  the 
radar  reflectivity  and  mean  vertical  velocity  profiles  [Lhermitte  and  Atlas,  1963], 

The  profiles  in  Figs.  6  and  7  exhibit  the  same  increase  of  mean  Doppler 
velocity  associated  with  acceleration  of  particles  due  to  melting,  wrhich  was  oh 
served  with  the  centimeter  wave  radars.  There  is  also  a  local  increase  of  radar 
reflectivity  at  the  top  of  the  melting  layer  attributed  to  the  change  of  index  of 
refraction  from  ice  to  liquid  wrater.  However,  the  decrease  of  signal  intensity 
below  the  bright  band,  which  was  always  observed  in  stratiform  rain  conditions 
with  centimeter  wave  radars,  does  not  occur.  This  can  be  seen  in  more  detail  in 
Fig.  8  which  shows  a  vertically  expanded  presentation  of  the  profiles  shown  in 

Fig  .6. 
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The  onset  of  radar  reflectivity  increase  occurs  slightly  above  (less  than  30m) 
the  height  at  which  particles  start  to  accelerate  which  is  consistent  with  cen¬ 
timeter  wave  observations.  However,  as  mentioned  above,  the  radar  reflectivity 
at  94  GHz  does  not  decrease  in  the  more  actively  melting  region  below  as  was 
observed  with  centimeter  wave  radars.  Also,  the  reflectivity  structure  above  the 
melting  zone  appears  as  a  notch  in  the  reflectivity  profile.  Indeed,  moving  down 
from  3800  meters,  while  the  particles’  mean  Doppler  keeps  increasing,  certainly 
due  to  particle  growth  by  accretion,  the  reflectivity  decreases  to  a  minimum. 
The  reflectivity  minimum  thus  occurs  in  a  region  of  steady  growth,  indicated  by 
the  downward  velocity  increase  in  the  mean  Doppler  profile,  in  which,  within 
the  Rayleigh  assumption,  radar  reflectivity  should  continuously  increase  as  it  is 
systematically  observed  at  centimeter  waves.  The  discrepancy  may  be  explained 
considering  the  structure  of  the  Mie  function  at  94  GHz  shown  in  Fig.  3,  which 
exhibits  a  well  defined  backseat  ten  ng  cross  section  minimum  for  ice  spheres  of 
1.6  mm  diameter.  Considering  the  Mie  functions  in  Fig.  3.  this  may  imply  that 
the  particles  grow  roughly  from  a  1  mm  to  a  2  mm  equivalent  size  in  the  region 
just  above  the  melting  level. 

There  was  no  available  measurement  of  rain  intensity  at  the  radar  site  but, 
according  to  a  relationship  between  mean  Doppler  and  rain  intensity  proposed  by 
Lhermitte  [1990],  the  precipitation  intensity  was  estimated  to  be  a  few  mm  hr -1 
(mean  Doppler:  4  m  s_1).  The  probability  of  raindrops  sizes  greater  than  1  to 
2  mm  in  such  weak  precipitation  is  very  small.  Indeed,  the  ’dark  band'  above 
the  melting  level  is  always  observed  at  94  GHz  in  weak  precipitation  intensity 
and  disappears  only  in  moderate  to  heavy  rain  conditions.  Fig.  9  shows  the 
time  variation  of  vertical  profiles  in  t lie  form  of  reflectivity  contours  presented 
in  altitude-time  coordinates  for  a  time  interval  of  approximately  3  minutes.  Fig. 
10  shows  the  mean  Doppler  contours  during  the  same  period  of  time.  One  sees 
that  the  reflectivity  minimum  is  deepest  between  10:28:00  and  10:29:30.  The 
minimum  fills  as  the  precipitation  intensities  as  indicated  by  the  increase  of  the 
mean  Doppler  to  4.5  m  s-’  at  10:29:00. 

5.  Second  experiment 

It  was  recognized  that  the  Doppler  data  interpretation  would  be  drastically 
improved  by  observing  the  complete  Doppler  spectrum.  Therefore,  in  addition 
to  the  pulse  pair  processor  providing  signal  intensity  and  mean  Doppler  profiles 
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at  range  gates  spaced  by  400  ns.  a  new  unit  was  added  for  the  recording  of  the 
I  and  Q  Doppler  signals  at  a  one  range  gate  position.  As  mentioned  in  section 
2,  the  gate  could  be  systematically  moved  in  10  steps  of  200  ns  (30  m),  so  that 
continuous  profiles  could  be  obtained  at  selected  4S0  m  altitude  intervals  in  the 
precipitation.  The  unit  was  capable  of  recording  4096  S- bit  complex  (I  and  Q  ) 
samples.  The  FFT  processing  needed  for  spectrum  calculation  was  planned  to 
be  done  off-line  by  a  microcomputer  also  including  programs  for  the  analysis  and 
interpretation  of  the  data. 

The  94  GHz  radar  was  operated  inside  the  shelter  from  November  26  to 
December  16.  During  the  three  weeks  of  radar  operation,  various  clouds  and 
precipitation  were  observed,  but  stratiform  rain  observations  were  obtained  ( 1  to 
20  mm  hr -1  rain)  only  on  November  30  and  December  11,  19S7.  For  both  days, 
observations  at  vertical  incidence  were  conducted  jointly  with  the  GL  TPQ-11 
35  GHz  radar  and  the  94  GHz  Doppler  radar.  In  addition,  rain  intensity  and 
dropsize  distribution  measurements  were  continuously  made  by  a  raingage  and  a 
disdrometer  installed  in  close  proximity  to  the  94  GHz  and  35  GHz  radars. 

As  an  example  of  the  30  November  observations,  ten  successive  signal  inten¬ 
sity  and  mean  Doppler  vertical  profiles,  acquired  with  the  94  GHz  Doppler  radar 
5  s  apart,  are  shown  in  Fig.  11.  Signal  intensity  is  1  / R2  corrected  but  is  not 
corrected  for  absorption. 

The  signal  intensity  profiles  exhibit  a  remarkable  time  stability  and  a  con¬ 
stant  l.SdB /km  slope  which  is  attributed  to  signal  attenuation  (two-way)  when 
propagating  in  the  rain  region.  The  slope  does  not  vary  appreciably  through 
the  vertical  profile  below  the  melting,  which  indicates  no  significant  evolution  of 
the  dropsize  distribution  from  the  melting  level  to  the  ground.  This,  together 
with  the  time  stability  of  the  profiles,  indicates  steady  and  homogeneous  rain 
conditions  between  the  melting  level  and  the  earth's  surface. 

Such  data  acquired  in  steady  stratiform  rain  conditions  provide  an  opportu¬ 
nity  for  accurate  measurement  of  signal  attenuation  at  millimeter  wnv«  and  it* 
relation  to  rain  intensity.  The  most  reliable  rain  attenuation  rneasuiement  w;,., 
made  in  an  S  mm  lir~l  rain  and  yielded  a  3.8  dB  km~x  value  which  is  shown  by 
the  circled  cross  in  Fig.5.  One  secs  that  this  data  point  agrees  well  with  other 
experimental  measurements  made  by  Wallace  [19SS]  also  shown  by  the  dashed 
line  in  Fig.  5. 

Fig.  11  shows  that  the  signal  intensity  variation  from  profile  to  profile  at  con¬ 
stant  altitude  is  very  small  (approximately  less  than  2  dB).  Using  a  relationship. 
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determined  earlier  for  94  GHz  radars  (Lhermitte,  1990],  between  mean  Doppler 
from  raindrops  falling  at  their  terminal  velocity  and  rain  intensity,  the  maxi¬ 
mum  mean  Doppler  variation  allowed  by  the  observed  signal  intensity  variation 
is:  Dons'1.  This  does  not  agree  with  the  observed  mean  Doppler  variability 
at  constant  altitude  (almost  1  tn  s~!  peak-to-peak)  which  must  be  therefore  at¬ 
tributed  to  up-downdrafts  reaching  approximately  ±50  cm  s~l.  This  shows  that 
even  generally  stratiform  conditions  such  as  observed  here  may  be  associated 
with  a  certain  amount  of  convection  in  the  rain  below  the  melting  level. 

The  systematic  increase  of  vertical  velocity  with  increased  altitude  is  due  to 
the  change  of  air  density  and  can  be  used  to  determine  the  exponent  in  the 
relationship  proposed  by  [Foote  an  DuToit.  1909]: 

1(D)-  1  of  D  i  ( iic./fh  V  ( 4 ) 

The  variation  with  altitude  of  t).<  a>  -  rage  mean  velocity  fnnu  data  shown  in 
Fig.  11  yields  x  ~  0.45. 

Comparisons  were  made  between  i!i<  signal  intensity  profiles  observed  at  94 
GHz  and  35  GHz.  An  example  of  the  ivmlts  is  shown  by  t!;e  three  successive  ft1 
normalized  profiles  observed  T. :  i  i  i ;  ■  ;  i,  •  -.o  ,e  »i:t.e  h.'ei  Town  in  Fig. 12  (94 
GHz),  and  in  Fig. 13  (35  Gllzi.  The  profiles  are  not  drastically  different  except 
at  the  top  and  above  tie  melt:,;-  lay*  r,  where  the  local  increase  of  reflectivity 
associated  with,  melting  .  ,5'  d-  t>  nt  ( approximate  ly  10  dB  at  94  GHz  and 
12  dB  at  35  GHz.  attenuation  corrected),  and  the  fact  tliat  the  35GHz  data  do 
not  exhibit  the  weak  reflectivity  minimum  above  the  melting  zone  seen  in  the  94 
GHz  profiles. 

The  slopes  of  the  vertical  pi o tiles  at  both  35  GHZ  and  94  GHz  are  surprisingly 
similar.  Since  the  attenuation  coefficient  t!1  35  GHz  must  be  around  1.7</Z?Gn_1 
(one  way)  [Lhermitte.  d  !'  ‘in  tic  ol<-ei ve«l  S  mm  In  1  lainfall.  the  7  dB  km~l 
two-way)  slope  tit  35  GHz  must  be  in  ptirt  due  to  some  misalignment  of  the 
two  (transmit  and  receive)  TI’Q-il  antennae  Aligning  two  transmitting  and 
receiving  antenna  beam-  can  9  a  m.i  V’i  problem  which  was  carefully  considered 
for  tin  -  9  1  GHz  :  i'  !m  .  f.  -i  v  i  ’■  :  !■•  ■  .  e:-  ■>  me  c«.iii-  v  : :  !,:n  3e.>s  1 1 -  an  0.1'  . 

Comparing  91  GHz  and  c,  ;.i  ;;a<  t-  :  wave  imlar  pn >!i!es  may  constitute  a  pow- 
erful  tool  for  the  nioiiitoring  of  precipitate, n  particles  size  and  shape,  on  the  ha- 
sis  of  which  their  evolution,  wht  n  falling  through  the  melting  level,  from  an  ice 
crystal-snowflake  complex  geometry  to  the  ultimate  and  simple  raindrop  shape 
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can  be  investigated. 

A  sizeable  effort  was  applied  to  calculation  and  interpretation  of  Doppler 
spectra  from  I  and  Q  signals  recorded  in  the  manner  presented  above.  The 
spectra  were  computed  from  FFT  programs  implemented  on  a  microcomputer. 
For  each  spectrum,  a  4096  complex  sample  Fourier  transform  was  performed. 
This  was  followed  by  spectrum  smoothing  done  by  a  running  filter  with  S  cm  .s'  1 
bandwidth,  providing  a  smooth  spectrum  (approximately  20  degrees  of  freedom 
for  each  spectral  density  estimate),  but  still  preserving  an  adequate  frequency 
resolution  in  the  16  in  s_l  unambiguous  Doppler  velocity  interval  allowed  by  the 
pulse  repetition  rate. 

An  example  of  such  spectra  observed  at  vertical  incidence  on  30  November 
in  a  8  mm  hr~l  rain  is  shown  in  Fig.  14.  The  spectrum  reveals  a  clearly  defined 
spectrum  dip  associated  with  the  minimum  of  backscattering  cross  section  seen 
in  the  Mie  scattering  functions  in  Figs.  2  and  3.  Figs. 15  and  16  show  vertic  al 
profiles  of  Doppler  spectra  in  a  450m  altitude  interval,  displayed  in  the  form  of 
spectral  density  contours  in  velocity  vs  altitude  coordinates.  Although  there  are 
noticeable  changes  in  spectrum  structure  and  bandwidth,  the  spectral  minimum 
(and  its  position  on  the  velocity  axis)  remains  remarkably  well  defined  and  stable, 
which  indicates  only  a  small  contribution  from  tip  and  downdrafts.  Note  that 
in  the  bottom  of  Fig.  15  another  minimum  at  8.2  m  .s_1  and  another  minimum 
at  approximately  9  m  s-1  occur.  This  is  due  to  an  increase  of  the  rain  intensity 
there,  with  the  dropsize  distribution  now  extending  to  3  to  4  mm  sizes  (see  the 
'theoretical'  spectrum  in  Fig.  17  below). 

Doppler  spectra  observed  at  vertical  incidence  (such  as  shown  in  Fig.  14)  can 
be  simulated  by  assuming  a  M-P  [Marshall  and  Palmer.  1948]  dropsize  distribu¬ 
tion,  -V  =  N0cxp{  —  \  D)  and  a  terminal  velocity  vs  diameter,  V(D),  relationship 
for  raindrops  derived  from  the  G-K  [Gunn  and  Kinzer,  1949]  experimental  data. 
The  expression  [Lhermitte,  1988b,  1990]: 

\ '(D)  =  V;[l  -  exp(  — 6.S  D2  -  4.8S  D)\  (5) 

is  an  excellent  fit  to  the  G-K  data.  In  this  equation,  V0  =  023(po/p: )° 45  with  pv 
and  pz  being  the  air  density  at  the  ground  and  at  the  altitude  £  at  which  V(D)  is 
expressed,  respectively.  Since  falling  raindrops  exhibit  some  random  variability 
of  their  fall  velocity  due  to  complex  interactions  between  the  drops  and  their 
environment  and  also  the  possible  occurrence  of  air  turbulence,  any  mean  drop 
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terminal  velocity  expressed  by  the  above  V(D)  expression  must  be  associated 
with  a  probability  function  simulating  the  drop  vertical  velocity  instability.  This 
effect  on  the  spectra  can  be  reproduced  by  submitting  the  predicted  spectrum 
derived  from  the  mean  V(D)  relationship  to  a  running  filter.  This  is  equiva¬ 
lent  to  convolving  that  spectrum  with  a  (Gaussian  shape  for  instance)  spectrum 
simulating  the  random  velocity  effect. 

Fig.  17  shows  a  comparison  between  an  observed  spectrum  and  a  theoretical 
spectrum  based  on  a  M-P  dropsize  distribution  with  A  adjusted  for  best  fit  of 
the  two  maxima.  Also  shown  is  a  correlation  between  observed  and  calculated 
spectra  as  a  function  of  a  velocity  lag  introduced  to  simulate  the  influence  of 
mean  air  vertical  velocity.  Only  the  spectral  components  between  3  and  7.5 
m  s-1  are  used  to  calculate  the  correlation.  One  sees  that,  the  correlation  is 
maximum  near  a  zero  velocity  lag  which  indicates  a  very  small  contribution  from 
air  velocity  (smaller  than  5  cm  ).  However,  the  theoretical  spectrum  exhibits 
a  deeper  minimum.  This  is  attributed  to  the  effect  of  the  random  components  in 
the  raindrop’s  fall  speed  mentioned  above,  which  smears  the  observed  spectrum. 

As  mentioned  above,  convolving  the  theoretical  spectrum  with  a  spectral 
function  simulating  the  random  velocity  component  corrects  the  problem.  This 
is  included  in  Fig. IS  which  now  shows  an  excellent  lit  except  for  the  discrepancy 
in  the  high  velocity  (large  raindrop’s  diameter)  part  of  the  spectrum  which  is  at¬ 
tributed  to  a  departure  fiom  the  exponential  slope  implied  in  the  M-P  dropsize 
distribution.  More  refined  spectral  models  such  as  the  Gamma  function  or  log¬ 
normal  spectra  [Willis,  19SS;  Feingold  and  Levin,  19SG]  may  be  used  to  provide 
a  spectral  fit  closer  than  that  allowed  by  the  M-P  spectrum. 

Fig.  19  shows  another  example  of  a  comparison  between  theoretical  and  ob¬ 
served  spectra.  A  random  velocity  correction  lias  been  applied  and  the  minimum 
depth  in  the  two  functions  is  in  good  agreement.  Howevei  the  maximum  corre¬ 
lation  is  now  obtained  for  -Idem  s-1  velocity  lag.  revealing  the  occurrence  of  an 
updraft  of  this  magnitn.de.  Note  that  measuring  the  position  of  the  minimum 
on  the  velocity  axis  is  an  adequate  estimate  of  air  velocity  and  appear?  to  be  a 
better  substitute  for  the  lnrtiind  of  spectrum  lower  bound  proposed  by  Battan 
[1963]. 

The  observed  variability  of  spectral  density  estimates  within  a  Doppler  spec¬ 
trum  observed  at  vertical  incidence  in  rain  is  greater  than  the  statistical  vari¬ 
ability  expected  from  the  degree  of  smoothing  applied  to  the  data.  It  is  thus 
interesting  to  note  that  a  substantial  volume  of  rain  (  more  that  3000  nr)  ob¬ 
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served  by  the  Doppler  radar  during  a  relatively  short  dwell  time  (a  few  seconds), 
still  yields  a  Doppler  spectrum  with  significant  small  scale  ’natural’  variability 
when  expressed  in  0. 1-0.2  mm  diameter  intervals. 

Dropsize  spectrum  measurements  can  be  derived  from  94  GHz  Doppler  radar 
observations  of  rain  at  vertical  incidence  [Lhermitte,  19SSb;  Lhermitte,19S9], 
Considering  the  above  discussion,  the  following  approach  is  recommended: 

1.  Define  a  M-P  dropsize  distribution  with  A  adjusted  for  best  fit  of  the  two 
maximum  amplitudes. 

2.  Perform  a  velocity  lagged  correlation  to  calculate  the  velocity  shift  at¬ 
tributed  to  air  velocity  (detecting  the  spectrum  minimum  velocity  position 
may  be  sufficient). 

3.  Convolve  the  calculated  spectrum  with  a  gaussian  spectrum  simulating  the 
random  velocity  component,  adjusting  the  bandwidth  for  best  correlation 
between  the  calculated  spectrum  resulting  from  that  operation  and  the 
observed  spectrum  shifted  by  the  velocity  lag  determined  in  2. 

4.  Convolve  the  Mie  backscattering  cross  section  vs  raindrops’  velocity,  <r(V), 
function  with  a  Gaussian  spectrum  having  the  bandwidth  determined  in  3. 

5.  Divide  the  observed  spectrum  shifted  according  to  the  velocity  lag  deter¬ 
mined  in  2  by  the  cr(V)  function  determined  in  4  to  obtain  the  dropsize 
distribution  expressed  as  a  function  V. 

6.  Calculate  the  dropsize  distribution  N(D)  =  N(V)dD/dV  using  the  (equa¬ 
tion  (5))  terminal  velocity-size  relationship. 

The  procedure  above  can  be  followed  automatically  by  a  computer  program 
which,  in  addition  to  the  dropsize  distribution  spectrum  provides  an  estimate 
of  vertical  air  motion  and  possibly  some  indication  of  air  turbulence  if  it  is  a 
significant  contribution  to  the  random  velocity  component.  An  example  of  the 
results  is  shown  in  Fig.  20  where  the  radar  derived  N(D)  distribution  is  compared 
to  that  obtained  with  a  disdrometer. 
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6.  Evolution  of  spectra  through  the  melting  level 

Doppler  spectra  were  observed  in  the  melting  layer  by  attempting  to  scan  the 
range  gate  at  which  I  and  Q  signals  were  recorded.  An  example  of  the  results  is 
shown  in  Fig.  21.  The  altitude  of  the  melting  level  had  to  be  estimated  at  the 
time  of  the  observations  by  interpreting  the  signal  seen  on  an  oscilloscope  so  that 
there  was  some  uncertainty  in  locating  its  exact  altitude.  This  is  the  reason  why 
most  of  the  vertical  profile  in  Figs.  21  appears  to  be  in  the  snow  region  above 
the  melting  level.  The  melting  starts  at  180  m  above  the  lowest  altitude  (420m) 
covered  by  the  range  gate  scanning  and  thus  does  not  cover  the  total  melting 
level  depth  (250m  ),  but  is  sufficient  to  determine  the  Doppler  spectrum  evolution 
through  it.  The  data  show  a  steady  evolution  of  the  spectra  with  greater  vertical 
velocity  variation  for  the  larger  size  edge  of  the  spectra. 

7.  Conclusion 

This  report  presents  new  data  acquired  in  stratiform  rain  using  a  millimeter 
wave  Doppler  radar.  Most  of  the  effort  was  applied  to  conducting  the  field  ex¬ 
periments  which  have  generated  a  decent  amount  of  data  considering  the  limited 
opportunity  for  stratiform  precipitation  observations  allowed  by  weather  condi¬ 
tions.  The  raingage,  the  disdrometer  and  the  TPQ-11  radar  were  very  valuable 
additions  to  the  the  project. 

The  3-mm  wavelength  radar  offers  the  unique  opportunity  to  ’label’  parti¬ 
cle  size  from  the  well-defined  backseat tering  vs  diameter  oscillations  occurring 
well  within  the  raindrops’  size  spectrum  range.  Indeed,  in  stratiform  rain  con¬ 
ditions  the  occurrence  of  very  small  air  vertical  velocity  contributions  can  be 
observed,  and  when  the  contribution  from  vertical  air  velocity  is  determined,  re¬ 
liable  dropsize  distribution  estimates  can  be  made.  In  more  convective  weather, 
measurements  of  up-downdrafts  may  be  derived  from  the  data. 

More  extensive  observations  of  stratiform  rain  conditions  in  various  meteoro¬ 
logical  situations,  based  on  vertically  pointing  94  GHz  Doppler  radar  observations 
assisted  by  AVband  and  S-band  radius  and  ground-based  dropsize  distribution 
and  rain  intensity  measurements  are  needed.  This  would  ultimately  provide  a 
better  understanding  of  stratiform  precipitation  processes,  such  as  the  influence 
of  snowflakes  melting  on  the  dynamics  and  physics  of  the  rain  below. 
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Fig.2  Radar  cross  secion  of  spherical  raindrops  as  a  function  of  their  diameter  at  35  GHz 
and  94  GHz.  The  dashed  lines  indicate  the  Rayleigh  scattering  approximation. 
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Fig.4  94  GHz  absorption  coefficient  for  liquid  droplets  in  the  Rayleigh  region  as  a  func¬ 
tion  of  the  droplets’  temperature  and  normalized  to  liquid  water  content. 
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RAIN  INTENSITY  (mm/lir) 


Fig.  5  Attenuation  coefficient  as  a  function  of  rain  intensity  at  35GHz  and  94  GHz 
computed  from  Mie  extinction  cross  sections.  The  dashed  line  is  a  least  squares  fit  for  the 
experimental  data  quoted  by  Wallace.  The  circled  cross  is  a  measurement  derived  from  the 
data  in  this  report. 
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Fig.6  Vertical  profiles  of  signal  intensity  (corrected  for  R 2  and  attenuation  and  expressed 
in  dBZ  radar  reflectivity)  and  mean  Doppler  velocity  (not  corrected  for  air  density)  on  15 
May  1987  10:30  EDT. 


Fig. 7  Profiles  taken  5s  apart  of  signal  intensity  and  mean  Doppler  (not  corrected  for 
air  density)  on  15  May  1987,  around  10:40  EDT.  (signal  intensity  is  corrected  for  R2  and 
estimated  signal  attenuation.  At  the  earth's  surface.  0  dB  is  roughly  equivalent  to  -10  dBZ. 
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Fig.  9  Time-height  display  of  radar  reflectivity  contours  in  the  melting  zone.  Time  starts 
at  10:'2S:00  and  ends  at  10:30:30,  approximately.  The  solid  line  is  a  melting  layer 
reference  also  shown  in  Fig.  10  for  comparison  between  velocity  and  radar  reflectivity  data. 
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Fig.ll  Presentation  of  10  profiles,  derived  from  data  acquired  on  30  November  at  5s 
time  intervals,  of  signal  intensity  and  mean  Doppler.  Signal  intensity  is  corrected  for  R2 
but  not  for  attenuation  and  indicated  in  relative  units.  Note  the  well  defined  signal  intensity 
vs  height  slope  equivalent  to  SdB  km~l  two-way  attenuation.  Mean  Doppler  data  are  not 
corrected  for  air  density. 


Fig.  12  R2  corrected  vertical  profiles  of  signal  intensity  but  not  corrected  for  attenuation, 
obtained  5  s  apart  with  the  94  GHz  radar  on  30  November  1987.  Same  data  as  in  Fig.  11 
but  now  shown  for  observations  made  at  the  same  time  with  the  TPQ-11  and  shown  in  Fig. 
13.  The  melting  layer  is  slightly  below  2500m. 
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SPECTRAL  DENSITY  RELATIVE  UNITS 


Fig.  14  Doppler  spectrum  obtained  at  vertical  incidence  on  a  8mm/ir-1  rain  on  Novem¬ 
ber  30,  1987. 
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Fig.  17  Doppler  spectrum  observed  at  vertical  incidence  together  with  a  predicted  spec¬ 
trum  (smooth  function)  based  on  the  Mie  backscattering  function  and  a  M-P  dropsize 
distribution  with  A  adjusted  for  the  best  fit  of  the  two  maxima.  The  function  at  the  top  of 
the  figure  is  the  cross  correlation  between  the  two  spectra  as  a  function  of  the  velocity  lag 
indicated  by  the  top  scale. 
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Fig.  18  Example  of  observed  and  predicted  spectra  adjusting  A  and  introducing  the  mean 
velocity  and  random  velocity  corrections  (see  text). 
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